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Abstract: We show how the sum of the two average tau polarisations in the decay chain 
X2 ~~ * T\t — * TT Xi m minimal supersymmetry with conserved .R-parity can be measured 
at the LHC. This is accomplished by exploiting the polarisation dependence of the visible 
di-tau mass spectrum. Such a measurement provides information on the couplings of 
the involved SUSY particles and allows a more precise determination of the di-tau mass 
endpoint. If different tau decay modes can be distinguished, the polarisation and endpoint 
measurement can be improved even further. 
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1. Introduction 



If Supersymmetry (SUSY) is realised in Nature at the TeV scale, the Large Hadron Collider 
(LHC) has a great potential to copiously produce sparticles. After discovery of new physics, 
the determination of its properties will be an essential task. This objective requires to 
perform as many measurements as possible to pin down the model. The vast majority of 
previous studies of SUSY at the LHC focused on methods to extract information on the 
superpartner masses from LHC measurements [jl|, |2|. 

In this paper we present a technique 
to obtain information on the average po- 
larisation of r-leptons originating from 
the decay chain 




xl 



t, r, 



T n T f Xl- 



(1.1) 



Figure 1: signal process, (3, 7 G {R, L} 



The indices n and / denote near and far 
to distinguish the two taus according to 
their mother particle X2 an d t\, respectively. As discussed below, the r polarisation pro- 
vides valuable information about the couplings of Xi\ X2 an d t\ and is important for a 
precise reconstruction of the t\ mass. The properties of the t\ are of special importance 
for relic density calculations, if the t\ is the next-to-lightest supersymmetric particle and 
is close in mass to the lightest SUSY particle (LSP) ||. 



- 1 - 



In the following we assume imparity conservation. As a consequence, sparticles are 
produced pairwise, each finally decaying into a stable LSP, which escapes detection, im- 
peding the direct reconstruction of SUSY mass peaks. Furthermore, the existence of decay 
chain ( |1 . 1[ ) implies the mass hierarchy m^o > > m^o . We also assume that the consid- 
ered di-tau final state is not significantly contaminated by other SUSY decays (e. g. three- 
body decays of xg, X2 ~* ^2 T n ~* Tn T Jx°i or X3,4 ~> ^l T n ~* Tn T Jx°i)- 

We exploit the shape of the visible di-tau mass spectrum polarisation sensitive 

observable. Since the energy of the tau-lepton decay products in the lab system is sensitive 
to its polarisation P Q, the invariant mass distribution of the visible final state is also 
affected by the polarisations P n and Pf of the two taus r ra and Tf [||, ||. The invariant 
mass distribution of the undecayed taus, m TT , has a triangular shape with a sharp drop-off 
at a maximum mass value 



2 

m- 




(see dashed-dotted line in Figure ||). The endpoint m™ T ax of the undecayed di-tau mass 
(including the neutrinos) is equal to the endpoint m^? x of the visible di-tau mass, since 
there is always a certain probability that m TT acquires its maximum value and in both 
tau decays the complete tau energy is transferred to the visible decay products. As the 
reconstruction of SUSY mass peaks is impossible, the determination of endpoints like 
Eq. ( |1.2| ) offers a convenient way to obtain SUSY mass information. The considered tau 
final state is of special interest, as it is the only di-lepton final state providing information 
on rriri ■ 

Unfortunately, in case of tau-leptons the endpoint measurement is more involved than 
for stable 1 leptons, since the visible mass spectrum m TT slowly fades out and gets very flat 
at large masses due to the escaping neutrinos from tau decays. Previous studies have shown 
that this problem can be overcome by choosing a suitable endpoint sensitive observable [l| , 
e. g. the inflection point of the trailing edge of the mass spectrum. It turned out, that even 
for low integrated luminosities, the impact of tau polarisation on this observable can already 
exceed its statistical uncertainty Q. In this case, polarisation is the largest systematical 
uncertainty on the m!^ x measurement and any progress on the determination relies 
on a proper consideration of tau polarisation effects. 

On the other hand, taus are the only leptons providing information on their polarisa- 
tion. This makes taus a unique probe of new physics giving access to information which is 
inaccessible in other measurements. 

This paper is structured as follows: In Section ^ the dependence of the polarisation of 
near and far taus on the underlying SUSY parameters is summarised. Section || discusses 
the shapes of the visible di-tau mass spectra as function of the tau polarisation for different 
tau decay modes. Finally, Section ||| presents the results of a Monte Carlo study using a 
parametrised fast simulation of an LHC detector m. It is analysed how tau polarisation 



1 Hcrc stable means stable on the scale of the detector size. 
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information can be extracted from the visible spectra for two different experimental scenar- 
ios: With and without discrimination of the tau decay modes. As an illustration, Section 



4.3 discusses the possible implications of such a measurement for the stau mixing angle 



and the stau mass under the assumption, that the Xi an d X2 properties are already known 
from other measurements. 



2. Polarisation of the taus 

The polarisation of the considered taus is determined by the coupling constants between 
T/3, P £ {72, L}, ri, and x% These are given by igAj^, which are shown in Figure [j]. On 
tree level A^ 1/3 can be written in terms of the mixing parameters of the neutralinos and the 
stau |§]: 

TYl 1 

A] 1L = -— -iV; 3 sin0 f + + tan IV) cos f (2.1) 

v^m^/cos/? J V2 

A T jlR = mr N j3 cos 9 f - y/2Nji tan W sin Of , (2.2) 

v2mv^cos/3 

where j = 1(2) for the coupling involving Xi(X2)i ^ji are the entries of the neutralino mix- 
ing matrix in the B, W, and H u basis Q , tan j3 = (H u ) / (H^) is the ratio of the two Higgs 
vacuum expectation values and Of the stau mixing angle. Here, only the lighter t\ is con- 
sidered. Within the chosen parametrisation of r-mixing, | t\) = cos Of \ tl) + sin^f | tr), 
the unmixed case is for Of = tt/2 resulting in t\ = tr and T2 = tl- 



The first terms of Eqs. (2.1) and ( p.2j ) are the Yukawa couplings of the Hd con- 
tributions to the neutralino state, preferring left chiral taus, tl, for a tr dominated t\ 
(sin^f ~ 1, cos Of ~ 0) and vice versa. This is due to the chirality nipping nature of the 
Yukawa coupling of the HiGGSinos. The H u has no contribution, as it does not couple to 
the down-type tau. The other terms are the electroweak gauge couplings, which conserve 
chirality and therefore the coupling to right (left) chiral taus is proportional to the Tril) 
contribution sin Of (cos Of). As the P^-field only couples to left chiral particles, the coupling 
to right chiral taus has no sensitivity to the wino contribution Nj2- The coupling to the B 
contribution Nj\ is enhanced by a factor of two due to the double hypercharge carried by 
right chiral particles compared to left chiral particles. 

With these couplings, the average polarisation is given by 



P = ) (2 . 2 . ( 2 -3) 



A hR) +{ A J1L 

where j = 1(2) has to be chosen for the T/( n ). 

Due to the scalar nature of the T\ , no angular momentum information can be exchanged 
between near and far tau. Therefore the realised polarisation state of the far tau is inde- 
pendent of that from the near tau in one decay chain. This is leading to the chiralities of 
the two taus being uncorrelated. 
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3. Shapes of mass spectra 



In the following Section we discuss the impact of tau polarisation on the visible di-tau mass 
spectrum for different tau decay modes. Our line of argument starts from the following 
assumptions: 1. The di-tau mass spectrum for undecayed taus has the triangular shape 
given by phase space. 2. There are no spin correlations between near and far taus. 

3.1 The decay r — > irv T 

The simplest case is the decay r — > ttv t . Due to the fixed chirality of the is T , as well 
as momentum and angular momentum conservation, the pion has a preferred emission 
direction depending on the tau polarisation P. The angular distribution of the pion is 
described by j| 

i Hr_ i 

BR(r^vr^ T )-(l + Pcost?), (3.1) 



T T d cos($) 

where ■!? is the angle between the pion and tau momentum. The pion coming from right 
(left) chiral taus is preferentially emitted along (against) the tau momentum, leading to 
harder (softer) pions. This is reflected in the di-pion spectrum, if both taus decay via 
r — * iru T . The spectrum is given by (see Appendix [A]) 



N(jn n7r ) = Ami (P n ■ Pf) mm (mm + 4m 2 + 4) + 4 (l - m 2 ) 



+ 



+ {Pn + Pf) 



m 2 — 2 In m — 1 — In 2 m 



+ In m 



(3.2) 



where m = m^/m™* is the di-pion mass normalised to the endpoint. 

The polarisations P n and Pf of 
the near and far tau enter the spec- 
trum in two different ways: through 
the sum P n + Pf and the product 
P n ■ Pf. Spectra for different sums 
of the two polarisations are shown 
in Figure [2] together with bands indi- 
cating the possible variations due to 
the product terms. For the product 
terms only those polarisation values 
are considered which are compatible 
with the respective sum of the po- 
larisations. Obviously, the dominant 
polarisation effect can be described 
by the parameter P n + Pf, whereas 
the product P n Pf has only a subor- 
dinate contribution. 

For the cases P n + Pf = ±2, the 



all polarisations 




0.2 0.3 0.4 0.5 0.6 0.7 0, 



Figure 2: m T 

product is fully characterised by the sum, whereas for the case of P n + Pf = the varia- 
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tion due to P n ■ Pf is most striking as here the individual pairs of P n and Pf satisfying 
P n + Pf = can vary most, e. g. P n = Pf = or P n = ±1 and Pf = =pl. 

In principle, a measurement of the m^-spectrum allows the determination of the 
endpoint, as well as the sum and the product of the polarisations P n and Pf, although, an 



extraction of P n Pf is challenging, as the effects due to the product terms in Eq. I3J2J are 
small compared to the contributions due to P n + Pf. If the product could be measured in 
addition to the sum, both polarisations could be measured individually up to the two-fold 
ambiguity of exchanging P n ^ Pf- 

3.2 The decays r —* pv T and r — > a\v T 

For the decay of the tau via the p and a\ vector meson resonance, the situation is more 
complicated as the longitudinal and transverse polarisation states of the vector mesons 
have to be distinguished. For the vector meson case, the fragmentation function is rather 
complicated j|, making the calculation of the resulting spectra cumbersome. Anyhow, the 
qualitative behaviour of these decays with respect to the polarisation can be described by 
the angular distribution of the decay products j| 

' ' '' - \ BR(t -> vv T ) , 2m l A l-P cos 0) (3.3) 



T T d cos i? 2 T 777,2 _j_ 2 m 2 

1 Tr 1 7?7 2 

-BR(t^w t ) „ ,X , (l + Pcos^), (3.4) 



T T d cos $ 2 + 2ml 

where T (L) denotes the transverse (longitudinal) vector meson state, P is the tau polar- 
isation, and $ is the angle between the vector meson and the tau momentum. Whereas 
longitudinally polarised vector mesons have the same angular behaviour as pions, trans- 
versely polarised vector mesons behave the opposite way, so that the expected net effect 
depends on the relative branching ratios between transverse and longitudinal states. These 
are described by the prefactors of Eqs. (|3.3j ) and (|3.4j) . With the values m T = 1777 MeV, 
m ai = 1230 MeV, and m p = 776 MeV [1(|, this leads to approximately the same amount of 
transverse and longitudinal states for the a\ resonance and more longitudinal ones for the 
p. Therefore, the mass spectrum for the decay via a\ is independent of the polarisation, 
whereas right (left) chiral taus have harder (softer) fragmentation for the decay via a p 
meson (similar to 7Ti/ T -decays). 

Again the spectra can be described by the sum and the product of the polarisations 
as can be seen from the following argumentation. For a given mean polarisation P the 
probability for a left-chiral tau is 5(1 — P) and that for a right-chiral tau is J(l + P)- 
Therefore the resulting spectra can be described by 

N (P n , Pf) = \{l + P n )(l + P f ) [RR] + i(l - P n )(l - Pf) [LL] + 

+ Pn)(l ~ Pf) + (1 " Pn)(l + Pf)} [RL] , (3.5) 

if [RR] , [LL] , and [RL] are the spectra for two right-chiral taus, two left-chiral taus and 



the mixed case, respectively. Eq. (3J3) leads to the same structure as Eq. ( |3.2[ ), beside the 
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Figure 3: Visible mass spectra for different polarisations and different decay modes. In each plot 



both taus decayed in the same channel. The curves are the functions of Eq. 3.2 



different distributions being proportional to the polarisation terms P n + Pf and P n ■ Pf. 
This argument holds in general as long as the there are no spin correlations such that the 
probability for the realised combined polarisation state can be factorised for the two taus. 

The visible invariant mass of the two taus for different decay modes is shown in Figure || 
for the cases of P n = P f = 1 (RR), P n = P f = -1 (LL) and P n = ±1, P f = =f! (RL) as 



obtained from the generator. The curves are the theoretical shape of Eq. (3.2) for the 
decay into single pions. In each plot both taus decay via the same channel. As expected, 
the decay over the a\ resonance yields a spectrum being approximately independent of 
polarisation, whereas the decays into single pions and over the p resonance are affected in 
a similar manner. 

In summary, the spectra can be described to a good approximation by the parameters 
m max anc j sum of the polarisations, P n + Pf, and the impact of the polarisation depends 
on the decay mode. 



4. Monte Carlo studies 



In the following we present an analysis method which allows to extract m^? x and P n + Pf 
from the visible di-tau mass spectrum. The feasibility for such a measurement is il- 
lustrated for the mSUGRA bulk region scenario SU3 (Mi/ 2 = 300 GeV, M = 100 GeV, 
Aq = —300 GeV, tan /3=6 and sgn(/i) = +1, Q) using an integrated luminosity of 36 fb _1 of 



pp collisions at y/s = 14 TeV. The SUSY sample is generated with Herwig v. 6. 510 uM, 12 



13], where the polarisation dependent tau decay is performed by Tauola [14, |lq| . QCD, 
tt, W+Jets and Z+Jets background are generated with ALPGEN v.2.06 pf. To simulate 
the response of an LHC-type detector, a parametrised fast detector simulation is used 
as an example. 

Due to the simplifications of the fast detector simulation and various idealising as- 
sumptions on the performance of the detector and reconstruction algorithms made in the 
following, the obtained results should be considered as a demonstration of feasibility rather 
than a precise evaluation of the expected detector performance. A detailed experimental 
analysis is beyond the scope of this study. A study using full simulation of the ATLAS 
detector showed that the decay chain of eq. [Hi] can be selected with high purity [||. 
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Figure 4: Mass and polarisation effects on detector level after reconstruction. The spectra are 
normalised to unity. Whereas the polarisation effects show up most obvious nearby the maximum, 
different endpoints can be seen mainly in the high mass part of the distributions. This will be used 
to disentangle mass and polarisation effects. 



We extend the technique used in Q, such that it also includes polarisation effects. To 
achieve this, we choose two observables, one of which ideally being maximally sensitive 
to the endpoint, the other one maximally sensitive to P n + Pf. The relation between 
these observables and the endpoint and the sum of the polarisations is established through 
a calibration procedure. The calibration is performed using 25 SUSY reference samples, 
where every combination of m™ ax = 80, 90, 100, 110, 120 GeV and P„ + P f = -2, -1, 0, 1, 2 
is realised. For these samples, only m^o, rn^o, m^, P n and Pf are changed to obtain the 
corresponding endpoint and sum of the polarisations, whereas all other parameters are 
fixed to the nominal values of the SU3 benchmark point. 

For suppression of standard model background, typical SUSY cuts ^ 230 GeV, 
Fr,i.jct > 220 GeV, PT,2.jct > 50 GeV, PT,4Jet > 40 GeV) are applied. Due to the large 
amount of missing transverse energy and the large number of hard jets in SUSY events as 
well as the two opposite sign tau signature, contributions due to background are negligible 

Q- 

To reconstruct the di-tau mass spectrum, the invariant mass of all reconstructed op- 
positely charged taus in one event is built for all events passing the selection cuts. Contri- 
butions due to fake taus 2 and taus not originating from the signal process (combinatorial 
background), can be efficiently suppressed by subtracting the corresponding distribution 
for same sign taus. This is possible, since the combination of uncorrelated taus has the 
same probability for equally and oppositely charged taus. Therefore, the background of 
the opposite sign (OS) distribution due to fakes and wrongly combined taus has the same 
shape as that of the same sign (SS) taus, since background with correlated OS tau-lepton 
pairs typically has small cross sections (e. g. xfxi production). 



2 E. g. non-tau jets, being identified as a tau-lepton. 
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Figure 5: The observables max (left) and Oy w (right) in the parameter plane of 
Pn + Pf- The two-dimensional histograms show the expected value of the observables for each 
simulated reference point (endpoint, P n + Pf), the red meshes represent the fitted function of eq. 



4.1 and the blue lines correspond to equal values of the observables. 



4.1 Without distinguishing tau decay modes 

Figure || shows the OS-SS mass spectrum for various combinations of endpoint and polar- 
isation, if no distinction between different tau decay modes is made. It is apparent that 
different polarisation states are distinguishable, although the discrimination power is re- 
duced compared to the di-pion spectra of Figure This is due to the fact, that the decay 
r — > ttv t is affected most by polarisation effects. Moreover some discrimination power in 
the low mass part of the spectra is lost, as the tau reconstruction is less efficient for low 
pr-taus. 

Figure ^| also shows, that SUSY masses and tau polarisation have different impact on 
the spectrum. Whereas the polarisation does not affect the endpoint and mostly shifts the 
position of the maximum, different SUSY masses stretch the whole spectra and show up 
most strikingly in the high mass part of the distribution. This will be used to disentangle 
polarisation and mass effects. 

For all 25 generated reference spectra, a fit with a Gaussian 3 is performed by Minuit 



1 17] and from this the position of the maximum, O max , as well as the position, where the fit 
decreased to b/io of its maximum value at the high mass part, Ou 10 , is calculated. Since the 
low mass part of the spectrum is strongly affected by the effects of the tau reconstruction, 
which does not contain any information about SUSY, the lower fit range is chosen to be 
the first bin, exceeding half of the maximal bin content of the spectra. This guarantees 
a common prescription for fitting the spectra, which is necessary in order to allow better 
comparability. However, due to statistical fluctuations of the bin content, the chosen range 
can vary, which needs to be considered in the error calculation. 

The 25 generated reference samples represent a 5x5 grid in the parameter plane 
spanned by m™ and P n + Pf. For each grid point the chosen observables C max and 



3 Among many tested functions, a Gaussian turned out to be the best compromise between shape agree- 
ment and mathematical simplicity. 
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Oy w and their errors can be calculated from the fitted Gaussian parameters. Figure || 
shows the obtained values for the two observables C max (left plot) and Oy 10 (right plot) 
plotted against the endpoint, m™ T ax , and P n + Pf. 

To obtain a continuous mapping from uv^^ and P n -\- Pf to the observables Omax 

and 

Oy 10 , a two dimensional fit of the function 

O ht (P n + P f ,m™*) = po {m™ x ~ 80 GeV) + Pl (P n + P f + 2) + 

+ P2 (m™ ax - 80 GeV) (P n + P f + 2)+p 3 , (4.1) 

is performed to the grid points for each observable. The fitted functions for the two 
observables are also included in Figure |5| (shown in red) . Both functions increase for larger 
endpoints, since the whole spectrum gets stretched for shifted endpoints. The dependence 
on P n + Pf is quite different: Whereas the polarisation has a sizable influence on the 
position of the maximum, Oy 10 , which is close to the polarisation independent endpoint, 
is hardly affected (compare Figure |||). This difference is reflected in diverse behaviour of 
the blue contour lines in left and right plot in Figure 

The measured value of one of the observables, max and 0i/ lo , can then be transferred 
to corresponding and P n + Pf values using the calibration. Each of the two measured 

observables defines a contour line in Figure [|, which provides a set of m^ x and P n + Pf 
values, being compatible with the measurement. The intersection point of these two contour 
lines finally yields both, the endpoint position and the sum of the two polarisations. It 
is crucial for this method, that the two observables have different dependencies on the 
endpoint and the polarisation, giving two intersecting contour lines. 

To test the described method, the nominal SU3 spectrum is generated in order to 
determine the endpoint and the sum of the two polarisations. The two observables are ex- 
tracted as described above. Statistical uncertainties are estimated by Minuit. Systematic 
uncertainties due to the variation of the fit range are determined by calculating the ob- 
servables for different choices of the fit range. Both contributions are added quadratically, 
neglecting possible correlations. 

To find the combination of m^? x and P n + Pf, fitting best to the measured observables, 
0™ and O^ e as -, the function 



T 



X z (m™ ax , P n + Pf) = \ O bt - O mcas -j Cov" 1 [Oj \ O m - 0™j (4.2) 

is minimised, where O fit = (0^ ax , Oy 10 ) are the expected observables values for given 
m^ ax and P n + Pf, as obtained from Eq. ( [4.1] ) and Cov(O) is the corresponding covariance 
matrix. This matrix can be calculated from the covariance matrix of the parameters of 
the fitted Gaussian. This method ensures the correct treatment of the correlation between 
Cmax and Oi/ 10 - The uncertainties on m^_ ax and P n +Pf are given by the contours, satisfying 
Ax 2 = Xmin + o(CL), where a(CL) depends on the desired confidence level. 

The result can be seen in the left plot of Figure |6|, where contours for different con- 
fidence levels together with the true value of the investigated SU3 point are shown. The 
determination of the endpoint is rather precise, whereas the polarisation is less constrained 
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contributions due to uncertainties in the cali- 
y mode is made. 



by the measurement. The tilt of the contours reflects the anti-correlation of the endpoint 
and the sum of the polarisations. As both observables increase with growing m™ ax or 
P n + Pf, a larger m^? x is more compatible with a smaller P n + Pf and vice versa for a 
fixed observable. 



In the above analysis uncertainties on the parameters pi in Eq. 4.1 are not considered 



Due to the limited statistics in the calibration samples, the estimated values for pi in 



Eq. |4.l| have errors, involving an additional uncertainty, which may be reduced in principle. 
Their contributions to the uncertainty can be seen in the right plot of Figure ||. Here the 
central value lines obtained from the measured values of O max and Oy 10 are shown with 
corresponding uncertainty bands. These bands are determined by solving 

O fit (m^ ax , P n + P f ) = O meas - (4.3) 

for P n + Pf and calculating the derivatives 



u y (x) 



\uU dlHdPj 



where y = P n +Pf and x = m^ ax . The effects are small compared to the other uncertainties 
and are neglected. The intersection of the lines is the central value of the measurement, 
since here x 2 = 0- 

4.2 With distinguishing tau decay modes 

So far, the shown spectra include all tau decay modes. As shown in the Section |3|, 
the polarisation effects depend on the decay mode. The decay via the a\ vector meson 
e. g. has no sensitivity on the polarisation and therefore the most precise determination of 
the endpoint can be achieved with events, where both taus decay via the a\ resonance. In 
contrast to that, the decay r — > itv T is affected most and is therefore a perfect polarisation 
analyser. Exploiting these differences between decay modes can improve the precision of 
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the combined endpoint and polarisation measurement. In the following we investigate how 
much improvement can be gained with a more sophisticated analysis taking decay mode 
specific differences into account. For this study we assume an idealised diagonal efficiency 
matrix for tau reconstruction: 



'akcs 



truth 



/BR(r -> X x y 



\BR(r X n ) 



(4.5) 



s 2 
i 1-5 
°- 1 



0.5: 
0^ 
-0.5; 

-1 
-1.5P 



39.35 % Confidence Level 
86.65 % Confidence Level 
98.98 % Confidence Level 
estimated value 
SU3 value 



where e = e(pT,f]) is the reconstruction efficiency averaged over all tau decay modes and 
iVr^x, are the number of taus decaying to X, L . Awakes is the number of fake tau-leptons 
from the parametrised tau-reconstruction algorithm. Fake taus are randomly assigned to a 
certain decay mode according to the respective branching ratio. The mean reconstruction 
efficiency for hadronically decaying taus with < \r/\ < 2.5 and 10 GeV < p™' < 100 GeV 
is 0.33. 

Now we analyse the di-tau mass 
spectrum separately for each tau de- 
cay category. In order to have rea- 
sonable statistics in each spectrum, 
splitting up all reconstructed di-tau- 
events into too many distributions is 
not viable. Hence, only two types of 
decays are differentiated: decays into 
single 7r, single K and p, which are 
affected by polarisation, and decays 
which are not affected (i. e. decays 
via a\ and all remaining decays). In 
this case one obtains three different 
cases: the fragmentation of both taus 
is polarisation independent, the fragmentation of both taus is polarisation dependent and 
the mixed case. Since for each of the three spectra C max and Oy 10 are measured, we end 
up with six observables in this case, each providing a contour line in (m^ x ,P n + Pf) space. 
The observables are determined in the same manner as for the case without distinguishing 
the tau decay modes. 

The result can be seen in Figure 0, showing, that the measurement can be significantly 
improved provided that information about the decay mode is available. Both the endpoint 
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Figure 7: Result with idealised tau decay mode dis- 
crimination 



and the polarisation measurement are affected. For these results the x 2 function of Eq. 4.2 
has been minimised using the six observables. 



4.3 Interpretation in terms of Of and rrif 



The measured parameters m^? x and P n + Pt depend on many fundamental parameters of 
the underlying SUSY model (Eq. |1.2| and 2.1-2^). Therefore, these quantities cannot be 
extracted in a direct way from the described measurement. Under the assumption, that the 
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Figure 8: The functions (P„ + P/)(0f ) and m™ 1 S lx (mf 1 ) with the measured values (red). 



neutralino sector (masses and mixing) is already known from other measurements, the sum 
of the polarisations is only a function of the stau mixing angle and the endpoint is com- 
pletely specified by the stau mass. In this case, the probability density /(m^ ax , (P n + Pf)) 
of Figure ||] can be translated into a probability density g(mf, Of) via 

g(m f , Of) = f (< n r ax (m?) , [P n + Pf] (0 f ) ) 

For that, the uncertainty contours of Figure |6] are approximated by ellipses making / a 
2-dimensional Gaussian probability density. 

For given neutralino masses the maximal possible endpoint of Eq. ( |1.2D is the difference 
of m.^o and m^o. Assuming the nominal neutralino masses for the SU3 scenario, this upper 
limit leads to the additional requirement mi^ ax < 101 GeV, which is more constraining 
than our measurement (see Figure ||). As this extra restriction is a special feature of 
the chosen benchmark point, we choose a different neutralino sector for our following 
analysis, where the kinematic limit on m^ ax is less constraining than the measurement. 
An example is the neutralino sector of the SU1 point (My 2 = 350 GeV, Mq = 70 GeV, 
Aq = GeV, tan/?=10 and sgnQu) = +1, [jij). For this scenario the functions m™ iB (mf) 
and (P n + Pf)(9f) are displayed in Figure || together with the expected precision of their 
measurement, shown as a Gaussian indicating their uncertainties. The stau mass must 
be in the range 139.5 GeV = m^o < rrif, < m^o = 262.9 GeV to allow the subsequent 
two-body decays, whereas in principle all mixing angles Of are possible. 

The transformed probability density can be seen in Figure ^. Since m^ ax (P n + Pf) 
is consistent with two values of rrif (Of), four regions are possible. The small preferred 
regions near Of ~ 1.5, 1.7 are due to the fact, that here, a small variation on Of causes a 
large variation in P n + Pf and therefore a large fraction of the possible values of P n +Pf > 
are projected on a small region of Of, making this values very likely. The long tail towards 



5m max 



d[P n + Pf 



dmf 



dOf 



(4.6) 
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3- 



higher (lower) values is due to the values P n + Pt < 0, which are distributed more equably 
on the #f-axis (compare left plot of Figure ||). 

The two possible are those be- 
ing consistent with the measured end- 
point (compare right plot of Figure |8|). 
The correlation of the parameters m™ x 
and P n + Pf can be seen in the bending 
of the four regions. 
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Figure 9: Of and mf determination 



5. Conclusions 

The polarisation dependence of the visi- 
ble di-tau invariant mass spectrum from 
the decay X2 ~~ * 



~o 

XiTT IS 



to a 

good approximation - parametrised by 
the sum P n + Pf of the near and far tau polarisation. It can be measured by exploiting 
shape variations for different polarisation states. Our method used two P n + Pf and m™ x 
sensitive observables to characterise the spectrum shape and establishes the relationship 
between the measured observables and P n + Pf and m^? x through a calibration procedure. 
From the polarisation measurement, information on the SUSY couplings can be derived, 



which might be valuable additional input for global SUSY parameter fits [18, 19, 

The precision on m^ x as obtained from a combined fit with P n + Pf can be improved 
with respect to an endpoint measurement neglecting polarisation effects. The results can 
even be improved provided that different tau decay modes can be distinguished sufficiently 
well. 
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A. Pion invariant mass distribution 



The invariant mass of the two undecayed taus can be written as 
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where •& is the angle between the momenta of the Tf and the t\ in the rest frame of the T\. 
As the T\ is a scalar particle, the probability density for fh TT a linearly increasing function. 

Within the relativistic limit, (3 T = 1, the probability F(x) of the fractional energy of 
the tau carried by the pion, x = E n /E T , is described by the fragmentation function 



Fi(xi) = [l + Pi(2 Xi -l)), 



(A.2) 



where i S {n, /} and P n j is the polarisation in terms of chirality. To obtain the resulting 
spectra of the two pions, the probability density H{z) for 



ml 



{Pir„ + Pn f ) 2 F 1Xn E nf (l- cos ip) E Wn E 7Tf 



mt 



[Pt 



+ Prf f E Tn E Tf (1-costf) E Tn E 



(A.3) 



is needed, where (p {•&) is the angle between the pions (taus). The fourth equality only 
holds in the collinear limit where ip = In the signal decay chain, the particle "between" 
the two taus is the scalar t, eliminating all spin correlations between the taus. Therefore, 
the probability distribution of z is the integral of F{x n ) ■ F(xf) over all x n and Xf with 
Z-n ' z f ^ z - The probability density can be obtained via the derivative 



H{z) = A ^5 dx n dx f F n (x n )- F f (x f ) 



x n -x f<Z 
xn,x f <1 



= {P n + P f ) lnz + 2(l-z) +{P n -P } ) 4(«-l)-Inz(4z + l) -lnz. (A.4) 

The di-pion spectra are the probability densities, N(m n7T ), of the variable = yfzrh TT . 
The probability density h(yfz) can be obtained from Eq. ( |A.4j ) via h{yfz) = 2yfzH{z) and 
that for m TT is the linear increasing function 2fh TT . With that, the pion spectra A r (m 7r7r ) 
can be calculated as 



N(m n7T ) = Am n7T [ -J-H (^f-) 
J m TT V m rr ) 



dm T 



(A.5) 



leading to Eq. 3.2 
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